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Abstract 



In the probability representation of quantum mechanics, quantum states 
are represented by a classical probability distribution, the marginal distri- 
bution function (MDF), whose time dependence is governed by a classical 
evolution equation. We find and explicitly solve, for a wide class of Hamil- 
tonians, new equations for the Green's function of such an equation, the 
so-called classical propagator. We elucidate the connection of the classical 
propagator to the quantum propagator for the density matrix and to the 
Green's function of the Schrodinger equation. Within the new description 
of quantum mechanics we give a definition of coherence solely in terms 
of properties of the MDF and we test the new definition recovering well 
known results. As an application, the forced parametric oscillator is con- 
sidered . Its classical and quantum propagator are found, together with 
the MDF for coherent and Fock states. 
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INTRODUCTION 



Quantum states are usually described in terms of wave functions |l| (for pure states) 
or by means of the density matrix 0, || (for mixed states). Nonetheless, since the 
beginning of quantum mechanics there have been attempts of understanding the notion 
of quantum states in terms of a classical approach Due to the Heisenberg prin- 

ciple, it is not possible to introduce a joint distribution function for both coordinates 
and momenta, while these joint probability distributions are the main tool to describe 
physical states in classical statistical mechanics. It is such kind of problems which has 
brought to the introduction of the so called quasi-probability distribution functions, 
such as the Wigner function ||, the Husimi function and the Glauber-Sudarshan 
function ||, |J, later on unified into a one-parametric family |fL0|| . Despite their wide 
use in quantum theory and their fundamental role in clarifying the link among classi- 
cal and quantum aspects, these quasi-probability distributions cannot play the role of 
classical distributions since, either they allow for negative values (like the Wigner func- 
tion) or they do not describe distributions of measurable variables. A formulation of 
quantum mechanics which is very similar to the classical stochastic approach has been 
presented by Moyal []TTj] . But, the evolution equation suggested by Moyal was an equa- 
tion for a quasi-probability distribution function (the Wigner function) and not for the 
probability. 



In [|12| , |T3l1 it was recently suggested to consider quantum dynamics as a classical 
stochastic process described namely by a probability distribution: the so called marginal 
distribution function (which was discussed in a general context in |10|), associated to the 
position coordinate, X, taking values in an ensemble of reference frames in the phase 
space. Such a classical probability distribution is shown to completely describe quantum 
states [ |PI| , |il3| | . The approach of JI2] - was developed both for quadrature observables 
@ and for spin 

Within this approach the notion of "measuring a quantum state" provides the usual 
"optical tomography approach" ||2l|-f23|j and its extension called the "symplectic tomog- 
raphy" formalism |p^,[15| . Both allow for an explicit link between the MDF and the 
Wigner function or the density matrix, in other representations. In this way, starting 
from the evolution equation for the density matrix, an evolution equation of the Fokker- 
Plank type for the marginal distribution function is obtained [12,13]. Such an equation 
allows for an independent definition of the marginal distribution function. Thus it may 
be assumed as the starting point for an alternative but equivalent formulation of quan- 
tum mechanics in what we call the probability representation, or classical-like description 
of quantum mechanics. In such a scheme, it plays the same role that the Schrodinger 
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equation plays in the usual approach to quantum mechanics. Since the MDF may be 
interpreted as a classical probability distribution, the Green's function connected to its 
evolution equation is called the classical propagator. In the classical propagator 



for a wide class of Hamiltonians and its relation to the quantum propagator for the 
density matrix are found. This establishes an important bridge among the probability 
representation of quantum mechanics and other formulations such as the path integral 
approach. 

The present work deals with two different problems. On one side we extend the 
emerging new description of quantum mechanics; on the other, we verify the formalism 
by applying it to a concrete, non-trivial case of physical interest, which is the forced 
parametric oscillator. A short report on the new results which we are going to present is 
already contained in |M] . From the point of view of the general formalism, we find new 



equations which connect the classical propagator of a quantum system to its integrals 
of motion. Also, we give a new definition of coherence, solely in terms of properties of 
the MDF, which relates the coherent marginal distribution function to invariants of the 
quantum system. In this regard let us note that the general approach to time-dependent 
invariants in quantum mechanics and their relation to the wave function were elucidated 



by Lewis and Riesenfeld in [25], while the connection among integrals of the motion and 



the quantum propagator was found in [g^,^]. In [f2HJ the relation of time-dependent 
invariants to the Schwinger action principle was established, and in the relation 
with the Nother theorem was discussed. 

Let us come to the forced parametric oscillator. This is a phenomenologically inter- 
esting model as it yields a good description of physical systems, like for example ions in 
a Paul trap. For such a system the symplectic tomography was already discussed in |18| 
while the endoscopy scheme for measuring states was suggested in fl3~0fl . The trapped 
ion may be in different nonclassical states like nonlinear coherent states ||31||. A general 



analysis of such states may be found in [32|. Here we try a complete description of the 



model in the framework of the new scheme. 

In section 1 we review the formulation of quantum mechanics in the probability repre- 
sentation. We introduce the marginal distribution function, together with its evolution 
equation. Then, we discuss the relation among the classical and quantum propagators 
and we specialize to the case of quadratic Hamiltonians. In section 2 we find two new 
equations for the classical propagator which is shown to be eigenfunction of a certain 
time-dependent invariant. We solve these equations for quadratic Hamiltonians and ob- 
tain the classical propagator as a function of time-dependent invariants of the system. 
This is an important step towards the characterization of the quantum system, as it 
yields both the quantum propagator for the density matrix (which we use to test the 
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scheme by comparing our results to well known results achieved by usual methods) and 
the time-dependence of the MDF, once it is known at t — 0. In section 3 we give a char- 
acterization of coherence directly in terms of the MDF, and we test the new approach by 
explicitly finding the coherent marginal distribution for the forced parametric oscillator 
both in the new framework and by means of more conventional methods. Finally we 
find the MDF for Fock states and we show a few significant plots. 



I. THE PROBABILITY REPRESENTATION OF QUANTUM MECHANICS 



In ref. |14]] an operator X is discussed as a generic linear combination of position and 
momentum operators 



X = pq + up 



(1.1) 



where p, u, are real parameters and X is hermitian, hence observable. The physical 
meaning of p, v is that they describe an ensemble of rotated and scaled reference frames, 
in classical phase space, in which the position X may be measured. In the above men- 
tioned paper it is shown that the quantum state of a system is completely determined 
if the classical probability distribution, w(X, p,u), for the variable X, is given in an 
ensemble of reference frames in the classical phase space. Such a function, also known 
as the marginal distribution function, belongs to a broad class of distributions which are 



determined as the Fourier transform of a characteristic function |H| . For the particular 
case of the variable (|TT|), considered in [12]- [15|, the scheme of [10| gives 



w(X,p,v) = — 
In 



dke~ ikX < e ikX > 



(1.2) 



where < A >= Tr(pA), and p is the density operator. In (10| it was shown that, 
whenever X is an observable, w(X,p,u) is indeed a probability distribution, as it is 
positive definite and satisfies the normalization condition 



J w(X, p, v)dX 



(1.3) 



The definition of the MDF allows us to express it in terms of the density matrix 

1 



w(X, p, v) 



2ttu 



p(Z, Z') exp 




dZdZ' . 



(1.4) 



Recalling the relation among the Wigner function and the density matrix, ( |I4| ) may be 
rewritten as a relation among w and the Wigner function , 
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w(X,p,u) = J exp[— ik(X — pq + vp)]W(q,p) 



dk dp dq 
(2vr) 2 



(1.5) 



Although the general class of distribution functions of the kind (|I.2|) was introduced, as a 
function of the density matrix, already by Cahill and Glauber in |10] , they didn't analyze 
the possibility of a new approach to quantum mechanics, in terms of such distribution 
functions, mainly because the invertibility of ( P~^ ) was not investigated. An important 
step in this direction is represented by f3]J where Vogel and Risken have shown that for 
a particular choice of the parameters p and v (the homodyne quadrature) the marginal 
distribution determines completely the Wigner function via Radon transform, namely 



they prove that (|I.5|) may be inverted for the Wigner function. In the same spirit, it was 



shown in |12| , |13| that the relation among w and the density matrix can be inverted for 
p yielding 



p{X, X') = -L / w(Y, p,X- X') exp 



Y -p 



X + X r 



dpdY 



(1.6) 



and, because of such a relation, the marginal distribution function satisfies an evolution 
equation 



d t w + Ow = 



(1.7) 



where O is a finite or infinite operator polynomial in q, p (also depending on p, u), 
determined by the Hamiltonian. Hence, we make our previous statements more precise, 
by saying that the MDF (the classical probability associated to the random variable X) 
contains the same information on a quantum system as the density matrix. The MDF 
may be defined through (|Oj ) independently from the density matrix, hence it represents 
the starting point for an alternative (but equivalent) approach to quantum mechanics, 
while ( |L~7| ) may be thought of as the analogue of the Schrodinger equation. Eq. ( |L~TD 
can be formally integrated to give 

w(X, p, v, t) = J n(A, /i, u } 0, X\ //, z/, t)w(X', p! ', z/', 0)dX'dpdu' , (1.8) 



where n(A, p, v, 0, X', p', z/, t) is the Green's function for the evolution equation (|L7|) . 
This is what we call the classical propagator. It can be interpreted as the classical 
transition probability density from an initial position X' in the ensemble of reference 



frames of the classical phase space, to the position X [13 



Let us now elucidate the connection of the classical propagator with the quantum 
propagator (Green function) for the density matrix p(X, X', t). Details may be found in 
I6| , [T9|| . For a pure state with wave function \I/(X, t), we have 
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p(X,X',t) = *(X,t)*\X',t). 



(1.9) 



Since the wave function at time t is connected to the one at initial time by the Green's 
function of the Schrodinger equation G(X,X',t) 

V(X,t) = [ V(X',0)G(X,X',t)dX' (1.10) 

we have for the density matrix 

p(X, X', t) = J K(X, X', Y, Y', t)p(Y, Y', 0)dYdY' (I.ll) 

with 

K(X,X',Y,Y',t) =G{X,Y,t)G*(X',Y',t). (1.12) 



The function K(X, X', Y, Y', t) is what is called the quantum propagator for the density 
matrix. Using the i 

K(X, X',Z, Z',t) 



matrix. Using the relation between the density matrix and the MDF ( |I.6|) we finally find 

1 



X + X r 



(2tt) 2 

x f — exp < % \Y — a 

J v' \ \ 2 J v 1 

x n(F, n, X - X', 0, Y', //, v\ t)dpdp'dYdY'dv' . 



7/2 



2v' 



(1.13) 



Then, once the classical propagator is known, the quantum propagator for the density 
matrix can be found. In the next sections, after finding the classical propagator, we will 
give an explicit example of this kind of calculation, for the driven parametric oscillator. 
We will also compare our results to analogous calculations obtained with the method of 
path integrals [ p3| . 



II. THE CLASSICAL PROPAGATOR 



In this section we address the problem of finding the classical propagator, in the 
framework of the time-dependent invariants method. We consider Hamiltonians of the 
form 

H = P - + V{q) (II.l) 

where V(q) is a generic potential energy. In this case the evolution equation for the 
MDF, Eq. (|T7|), takes the form [Q,0,|ll 
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d 

w — Utt w ~ 1 
dv 



V 



' -1 d 
d/dx dfi 




-Id v d x 
l_ i 

d/dx On 2<9x, 



w 



0. 



(II.2) 



Restoring the proper units, the Planck constant h will appear in Eq. ( [11. 2j ) so that the 
equation, even if classical-like, gives a quantum description of the system evolution, so 
replacing the Schrodinger equation in our scheme. It can be shown jl9| that in classical 
statistical mechanics the distribution w(X, /i, v) may be also introduced and the classical 
Boltzmann equation can be rewritten for this distribution. Then the classical limit of 



( [O ) is the Boltzmann equation. The classical propagator obeys an evolution equation 
which follows from ( |II.2| ) 



dU d „ 



V 



' -1 d 
d/dx d[i 




5{t 2 - t 1 )5(X' 



— Id v d 
l_ i 

d/dx dfi 2dx ; 
X)8{n' - fi)6(u' - 



n 



v) 



(II.3) 



with initial condition 



n(X, /i, v, 0, X', //, i/, 0) = 5{X' - X)5(fi' - fx)5(p' - v) 



(II.4) 



To be definite, let us consider the Hamiltonian describing the forced parametric oscil- 
lator. This is of the form described by ( |11.1| ), and it includes many other quadratic 
Hamiltonians as limiting cases. The potential, V(q), is given by 



v{q) = -^q 



fflq, 



(11.5) 



that is 



(II.6) 



We leave u free to be real or imaginary, to allow the description of repulsive oscillators as 
well. The evolution equation for the marginal distribution w, is obtained by the general 
expression (|TL2|) , replacing the potential by (|L5[) . We have 



w 



d ( 9 , . d . d \ 







(II.7) 



The equation for the propagator may be obtained analogously. 

The solution to ( [II. 3D was previously argued to be expressible in terms of the integrals 
of the motion of the system We show here a derivation of this result, which is 

essentially an extension of previous techniques p5|- |27l , p9| , |28| . In the mentioned papers 
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the time-dependent invariants of a given system were shown to be in connection with the 
wave function and with the Green function of the Schrodinger equation. The integrals 
of the motion, I(t) are defined by the equation 



^I(t)+i[H,I{t)] = 



(11.8) 



where H is the Hamiltonian of the system. We may think of time- dependent invariants 
as the evolution in time of the initial coordinates and momenta, q and p. Then we 
express them as 



J(t) = A(t)Q + A(t) 



(II.9) 



with Q = (p, q) and A, A to be determined. For the oscillator any integral of the motion 
may be expressed as a function of the two operators A(t) = UaU~ l and A^(t) = Ua^U^ 1 , 
where U is the evolution operator obtained by the Schrodinger equation. Of course, they 
are integrals of the motion in their turn. We have 



A{t) 
A\t) 



V2 



[e(t)p-e(t)q)+P(t) 



with e satisfying 



e + uj 2 {t)( 



ee — e e 



2i 



and initial conditions e(0) = 1, e(0) = i. It can be checked that 

[A, At] = 1 



(11.10) 
(11.11) 



(11.12) 
(11.13) 



(11.14) 



The function (3(t) is determined by consistency with ( |11.8| ) to be 

P(t) = -Asf *! <t')f(t'). 



The integral of the motion 



I(t) 



'A n A 12 
A21 A22 



p 



a/ 

A 2 



;n.i5) 



;n.i6) 



can be given in terms of A and A< as 
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_ A - A* 

OUT, 



The matrix Lambda and the vector A are then determined by comparing (|II.16|) with 
(020) • We obtain 

and 



In [p6| , p7| it was shown that the Green function G(q, q', t) is a solution of the system 

I q G(q,q',t) = q'G(q,q , ,t) (11.20) 
I p G(q,q',t) = -p'G(q,q',t), (11.21) 

where the first equation means that the Green's function is an eigenfunction of the 
invariant I g at each value of t, with eigenvalue the initial position, q', of the system. 
These results were derived by applying the evolution operator U (t) to the identities 

qd(q-q') = q'S(q-q') 

pS(q-q') = -p'S(q-q') , (11.22) 

where G(q,q',0) = 6(q — q'). Equations ( [II. 2Q ) and ( |II.21| ) may be trivially generalized 
to equations for the quantum propagator of the density matrix K(X, X', Y, Y', t). Thus 
in principle we can find analogous relations for the classical propagator II, inverting 
the relation ( |1.13| ). This procedure, although mathematically well posed is in practice 
difficult to pursue. We will use instead a formal procedure which goes along the same 
lines of the derivation for the Green's function. 



Assuming the existence of an evolution operator, U, for the equation ([1L7|) , and 



recalling that n(X, //, u, 0, X', //, u', 0) = 5{X - X')8(/i - yf)8{y - we find 

I q U(X, ii, v, 0, X', n', z/, t) = q'Il{X, n, v, 0, X', n' , i/, t) (11.23) 
I P U(X, ii, v, 0, X', it', v', t) = -p'U{X, ii, v, 0, X', ii', v', t) . (11.24) 

Let us explain the notation. We mean with <f and p' the operators which represent the 
action which is induced on w (X, /i, v, t) when acting on the density matrix with q' and p', 
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due to (0). Analogously I q and I p represent the action which is induced on w(X, p, u, t) 
when acting on the density matrix with I q and I p . The operators I q and I p are formally 
connected to q and p through the evolution operator U. They are given by 



q 



P' = - 



d \ 1 d v' d 



dX> J dp,' 2 dX' 







-i 



d .p! d 



dX' / dv> 2 dX 



; £ _ e .)__ (e _ £ .)_ 
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-1 



+ 



' y/2 

'J^Y 1 P -(3* 

dx + ^/T 



The solution to these equations can be verified to be 

U(X, p, u, 0, X', p', v\ t) = 5{X -X' + A/A^A^A"' - AA" 1 ) (11.25) 



where J\f, Af' are vectors, Af = (is, p), Af' = (zV ', p'). 

Substituting QII.18I) and ( |II.19| ) into (|II.25|) we get the classical propagator to be 



n( x, „, „, o, a", s,s, t )=i (x - r + »w+in+m±^ 



;n.26) 



/ „(e* - e) + p(e* - e) \ ( , v(e* + e) + + e) 
X 1 1/ — z o I // 



Now, we can replace this expression into ( |1.13| ) and we obtain the quantum propagator 
for the density matrix of the driven parametric oscillator. Of course, the integral cannot 
be performed unless we assign the explicit dependence of the frequency, u(t), so that we 
can solve ( |II.12| ) for e(t). Here we assume uj = constant. We have then 



e(t) = e u , 



(11.27) 



and the classical propagator ( |11.26| ) assumes the form 

n(A , p, u, 0, X , p , u , t) = d I X - X H -j= I 

x 5 iy' — (z/cost + /isint)) x 5 (p' + (z/sint — pcost)) . (11.28) 



Substituting into (|I.13|) , this yields the quantum propagator 



K(X, X', Z, Z\ t) = / I exp (i f y - - i——j—Y' + i^±^V 



(2tt) 2 7 i/ 



2i/ 
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x S [Y -Y 



, 0e- a ([i - i(X - X')) + (3*e lt (n + i(X - X')) 
x 8{v' - ((X - X')cost + fisint)) x 5 (// + ((X - X') sint - /icost)) d/idfi'dYdY'dv' . (11.29) 



Performing the integration we get 



K(X, X',Z, Z',t) 



27rsint 6XP 1 v^2 



(3e- u -i(X - X') + 



Z-Z' (X-X')cost 



x exp 



—i 

75 



P*e lt i(X - X) + 



sin t sin t 

Z-Z' (X-X')cost 



.(X + X>) 



x exp < —% 



x exp -(Z + Z') 



sin i sin i 

Z — Z' (X-X')cost' 



s'mt 



s'mt 



(11.30) 



, v . (Z-Z') cost (X-X')cos 2 t 
-(X - X')smt + - h 



sint 



sint 



where we have used 



dY' exp < iY' 



1 - 



:z - z') 



and 



jj sin t + (X — X') cos t 
1 



*(/(/*)) 



2tt5 1 - 



r 5(/x - flo), 



(Z - Z>) 



/zsint + (X — X') cost. 



l/'MI 

with /i s. t. /(/io) = 0. 

In view of interpreting the quantum propagator for the density matrix as the product 
of quantum propagators for the wave function, as in (|L12|) , we separate the primed and 
unprimed variables in (|II.30|) , obtaining 



K(X, Z, X', Z\ t) = - J— exp {^[(X 2 + Z 2 ) cost - 2XZ]\ 
27T suit 12 sin t J 



x exp 



—i 



sint 



[Z((3e- U + (3*e a ) - X((3 + (3*)} 



x exp |-^-^[(X /2 + Z /2 ) cost- 2X'Z']| 
[Z'(P*e a + f3e~ u ) - X'{p* + (3)] 



a/2 sin t 



;n.3i) 



This expression can be further simplified using the explicit form of the shift, (3, as given 
in (PT5Q . We finally get 



K(X, Z, X',Z',t) = — L- exp { -4— [(^ 2 + Z 2 ) cos t - 2XZ] 
27r sin t 1 2 sin t 

x exp ^2Z J f(lf) sin(t - t')dt' + 2X^ f(f) sint'dt' 
.[(X /2 + Z /2 )cost-2X'Z / ]| 



x exp 



2 sint 
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x exp S.2Z' J* f{t') sin(t - t')dt' + 2X' J* f{t') sint'dA . (11.32) 

Now we can read off the quantum propagator for the wave function. It is worth noting 
that it can be determined only up to a phase factor, independent on phase space, but 
possibly dependent on time, as can be argued from ( |1.12j ). We have 

1 



G(X,Z,t) = e 



iF(t) 



= cxp , 
V 2ir sin t 12 sin t 



x cxp 



( ' 


2Z f 


I 2 sin t 


Jo 



[(X 2 + Z 2 )cost-2XZ] 
et 



(11.33) 



where F(t) is an unknown, real function. This is the expected result, already obtained 
in the literature with other techniques (cfr. for example f3~3]| , where the quantum prop- 
agator is obtained with path integral methods). Thus, we have checked in a specific 
example that the probability representation gives equivalent predictions to the standard 
formulations of quantum mechanics. 

The quantum propagator for the simple harmonic oscillator (u =constant, /3 — 0) is 
easily obtained to be 

K(X, Z, X', Z', t) = — - — x exp ( —!— \(X 2 + Z 2 ) cos t - 2XZ] } 
2n sin t 12 sin t J 



x exp {-^-jl( X ' 2 + Z ' 2 ) cost - 2X'Z']\ 



(11.34) 



yielding 

G{X, Z,t) = e iF{t) - 



x 



exp(-^— [(X 2 + Z 2 )cost-2XZ] \ . (11.35) 



V 27T sin t 

The quantum propagator for the free motion, already calculated with the illustrated 



techniques in |T^,[T9[, may be recovered in the limit (3 = 0,lu = (which corresponds to 
sint — ► t, cost — > 1). We get 



K(X,X', Z, Z',t) = -Lexp {^[-(X - Z) 2 + (X' - Z') 2 ] 



[11.36) 



yielding 



G(X, Z,t) = e 



iF(t) . 



;n.37) 



12 



III. THE MARGINAL DISTRIBUTION FOR COHERENT STATES 



In this section we find the marginal distribution function for coherent states, using 
three techniques. First we give two derivations which are based on a reformulation of 
the notion of coherence in terms of the distribution function itself. These derivations 
are particularly relevant from a conceptual point of view giving another example of 
how common concepts of quantum mechanics are treated in the new approach without 
recursion to the wave function. Then we use the relation among the distribution function 
and the density matrix given by (|T6|). In this approach, the notion of coherence is 
defined in the conventional manner, through the wave functions, by requiring that they 
be solutions of the Schrodinger equation with initial condition 

atp a = aip a (III. 1) 

(with a, cfi usual annihilation and creation operators). As known the two equations (the 
Schrodinger equation and the initial condition) may be put together, yielding 

A(t)i) a = aip a . (III.2) 

where the operator A, defined in the previous section by ([Ll.llp , represents the time 
evolution of a. Once we find the solution of ( |111.2| ), we can write the density matrix 
associated to coherent states, and, consequently, the MDF, through (0|). 



A. Probability Representation Approach 

Let us describe the new approach, in the spirit of the probability representation of 
quantum mechanics. The first problem we are faced with, is to define the notion of coher- 
ence of a quantum state, independently from wave functions. We want coherent states 
to be defined as those states whose MDF satisfies the Fokker-Planck equation ( |11.2| ), 
with an initial condition which "translates" ( |IIL1|) into an equation for w. Applying the 
annihilation operator a to the density matrix, and using ( |111.1| ) we get 

ap(X, X') — ap(X, X') (III.3) 

(acting on the right with a) we get an equivalent equation, with a* instead than a). 
Now, recalling that a = (x + and using the expression of the density matrix in 
terms of the MDF, flOJ), Eq. (iHlll) yields 



1 (fi + iu d (d .d\ d - 1 } , v . , v . /jjj a \ 

^{ — M + [ty + t to)dX ]"a{X, l *,u) = a w a {X, l * t u), (III.4) 
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where we have used the correspondence J01 



Xp(X,X') 



1 _d_ 
2^dX 



-i 



9 > 








9 


i 9 



dX ) d/j + 2 U dX 



W(X, /!, V) 

w(X, /i, f) 



(III.5) 



Eq. (pl.4|) may be rewritten in the compact form 

aw a (X, fi, v) = aw a (X, /i, u), 



(III.6) 



where the operator 



1 ( u + iu d ( d d 

— < —~ — + 7T" + 



ax 



d_ 

dp) dX 



(III.7) 



translates the action of a on the density matrix into an action on w. Hence, we may 
reformulate the problem of finding the MDF for coherent states entirely in the language 
of classical probability, the initial condition fllll.lp being replaced by ( |ill.4[ ) . The solution 
of ( |ill.4| ) represents w a at t — 0. We get the solution at a generic value of t just as we 
would do for wave functions, namely, by applying the propagator for the corresponding 
evolution equation. For wave functions this is the Green function of the Schrodinger 
equation, while for the MDF it is the classical propagator which we found in the previous 
section (the Green function of the Fokker-Planck equation). To solve ([III.4j ) we make a 
Fourier transform of the marginal distribution function 



w a (X, /i, v) = I Wk(a, fi,is)e lkX dX 



(III.8) 



and we get 



ik 



fi + iu 1 



2 

This equation is of the form 



2 y + %, 



{-- 









y/2, 



OtWk 



(III.9) 



z 
2 



dz 



V2, 



(111.10) 



with y = kfi, z = kv. The solution of such an equation is known to be a Gaussian 

Wk ~ exp (cy 2 + dz 2 + hyz + ey + gz^j . (III. 11) 
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It turns out that, to determine the coefficients, we also need the complex conjugate of 

( prmp , 



' i d \ ( z d s 

2 y+ % + U"&. 



V2a*w k - 



(111.12) 



where we have used 



w t(y, z ) = w- k (y, z) = w k (-y, -z) 



(111.13) 



(remember that w a is real). Hence, we determine c,d,e,g,h by consistency to be 

(111.14) 



1 a - a 

c = a = — , e - 
4 



y/2 



Restoring the original notation we have 



.a* + a 
g = —% t=— , h — (J. 



V2 



Wk(ot, /i, v) = N exp < —k 



2 (/" 2 + ^ 



+ 



a — a 



V2 



kv — i 



a — a 



V2 



(111.15) 



where iV is a normalization factor. Taking the inverse Fourier transform of ([III. 15| ) we 
finally get 



w a (X, fi,u) = N 



7t(/J, 2 + V 2 



exp 



a+a* ,,\2 



V2 



/i 2 + V 2 



(111.16) 



The last step is to find the time dependence of the coherent marginal distribution func- 
tion. This can be achieved by substituting the solution at t — 0, ( |III.16|) , into Eq. (|L8|) , 
where the propagator is given by (|II.26|) . We have 



w a {X, fi, v,t) — N J II(X, n, is, 0, X', fjf, v', t)w a (X', //, is', 0)dX'dfi'du' 



N 5 \ X — X 



V2 J 



x 



7r(/i /2 + v r 



exp 



(X' + i2^V_2±fV) 2 ' 



dX'du'dv' 



fi' 2 + v' 2 



which, performing the integration, yield 

1 f 1 



w a (X,n, u,t) = N 



7r|ez/ + eji\ 



exp 



\eu + en\ 



X - 



7(e> + ve*) + 7*(e/i + ev) 



V2 



(111.17) 
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The normalization factor, N, is equal to 1, due to the normalization condition ( |1.3| ). 

It is worth noting that, whenever the time dependence is not too complicate, the 
previous calculation to get the time-dependent marginal distribution function, which 
consists of two steps, may be replaced by directly solving the time-dependent equation 

A(t)w a (X, /i, z/, t) = aw a (X, fi, v, t) (III. 18) 



which is obtained in the same way as ( |111.2| ) as a straightforward extension of the pro- 
cedure illustrated for the time independent case. The operator A(t), which we will find 
explicitly below is obtained from the action of A(t) on the density matrix when we ex- 
press p(X, X',t) as a function of w a . In our situation ( [111.18] ) turns out to be solvable 
using the same techniques which we used to solve the time independent analogue. This 
derivation is completely equivalent to the previous one. We choose to present both of 
them, because one requires the explicit use of the classical propagator, emphasizing the 
great content of information contained in it, while the other shows how simple calcula- 
tions can be in the framework of the probability representation of quantum mechanics. 

Since the operator A(t) given by ( [II. 10| ) may be rewritten in coordinate representation 
(of phase space) as 




ie(t)X)+(3(t), 



(111.19) 



we can use the correspondence (|III.5|) to write ( |IH.18|) as 

d \ ( d V 



1 J e/i + ev d . ( . d 
V2 [~^~ dX + \% 



du \dX 



w a (X, \x, v) = n/w a (X, fi, v) , (III. 20) 



where 7 = a — (3. As before, this can be easily solved performing a Fourier transform of 
the MDF. We get 



' i . d \ .fez d s 

,2 ty + % +i T- £ &. 



V2<yw k 



(111.21) 



with y = kfM, z = kv. The solution is a Gaussian of the form ( |111.11[ ). To determine the 
coefficients we proceed as before. We replace (|III.11|) into ( |III.21|) . We get a consistency 
equation for the coefficients, which doesn't determine them completely. Hence, we take 
the complex conjugate of (|IH.21| ), and we use the fact that w* k = w_k = Wk(—y, —z). 
We obtain 

— + 2ed - eh) y + (- - 2ec + eh) z + ee - ef = v^T 
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'?£_ _ 2e*d + e*h^j 2/ + ( Y + 2e*c - e*/i ) z - e*e + e*f = v 7 ^*, 



fIII.22) 



which yields 



C = ~X' d 



\e\ z _ 7e* + 7*e 
4 ' 6_ ' ; 



7e + 7 t 
z\/2 



. (111.23) 



Restoring the original notation we thus get 

f k 2 

u>k(a, fi, u, t) = N exp < — — /i 2 |e| 2 + z^ 2 |e| 2 + + (ee* + e*e)nv 



x 



4 



/ 7e +7 e . 7^ + 7 e 

-Z/C jj, 7= h V- 



s/2 



V2 



(111.24) 



Taking the Fourier transform of this expression, it is immediately verified that it is 
identical to ([III. 171 ), as expected. 



B. The Schrodinger Approach 

Let us come to the more conventional point of view. We first find the coherent wave 
functions and the corresponding density matrix, and then obtain the coherent marginal 
distribution function through the Eq. ( |1.4| ). The solution to Eq. (|111.2| ) is of the form 

i/) a (t)=C(t)exp(iY~ + ^(a-P)xJ • (111.25) 

We determine the overall factor, C(t), up to a phase factor, by imposing the wave 
function to be normalized to 1. We pose C(t) = D(t) exp(z</>), we have then 

1 = r \^ a \ 2 dx = D(t) 2 V^* exp (_ [( a - /3)e* + {a* - [3*)e} 2 \ , (111.26) 
J— oo 1 2ee* J 

from which we get 

(111.27) 

We now replace the solution we found in ([OJ), recalling that 

p Q (Z,Z')=MZW a (Z') . (111.28) 
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After some algebra we get 



w a (X,fi, v) 



1 



dZ exp 
dZ' exp 



exp 



-( 7 e*+ 7 *e) 2 



2ee 



x 



(111.29) 



where we have posed 7 = a — (3. The Z, Z' dependence is completely factorized and the 
two integrals can be put into Gaussian form so that 



w a (X,fi, v) 



zV ^ee* If + _ 

i /2 7 2 X 2 2^22 7 



x 



x exp 



x exp 



e 2 z/ 2 



X x 



z/ e 

2 7 * 2 X 2 2y/2i 7* 



2(^ + ^ 



+ 



(111.30) 



This rather complicated expression becomes simpler if we recognize the exponential 
factor as a square of three terms: 

2^ 



w a (X,/i,p) 



1 



vJnee* - + 



*| e 1 A* 1 2 



exp 



X 



v/2 



7^* (f + *) +7*e(f + £) 



zA/e*e - + 



> . 



(111.31) 



This expression can be seen to coincide with (|III.17|) , hence confirming the equivalence 
of the two approaches. By evaluating the mean value of X , 



< X >=< p > u+ < q > fj, 



(111.32) 



we find 



1 1 

<X>=-j= ( 7 e* + 7 *e) v+^= (7e* + 7*e) /i, 



(111.33) 



so that the marginal distribution function for coherent states ( |IH.31|) takes the simple 
form 



w a (X,n,v) 



1 



2na 2 x 



exp 



(x- < x >y 

2^1 



(111.34) 



where a\ is the variance of the variable X 



2 1 2 * I e 1 A* |2 
ffv = -1/ ee - H — 
x 2 'e z/ 1 



[111.35) 
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C. The Marginal Distribution Function for n-th Excited States 



From the expression of w a , ( |111.17| ), we may evaluate the MDF, w n , for the n-th 
excited state. We have 

1 { X 2 1 

w a (X,fi,u) = w ex P \ -\ a - P\ 2 - r- — ; — v? \ 

Vvr|ez/ + Lie\ 2 I \eu + Lie\ z \ 



x exp 
x exp 



2\€i' + e/x| 



|ef + efj,\ 



(111.36) 



which can be put into the form 

1 f , ,0 



w a (X,ri,v) 



X 



0r|r 
1 

\/7r|r 



exp < — \a\ 



\P\ 2 -V2—(r(3* +r*/3) 



2\r\ 2 | 



e -|a| 2 -r 2 e -i 2 +2ty e -t* 2 +2t*y 



(111.37) 



with 



r = ev + e/j, 
t = a- 



V2\ 



r 



Y = T2 



P*r + f3r* + ^2X 



Remembering the generating functions of Hermite polynomials: 

oo t n 



-t 2 +2tY 



E -,b»0O, 



n=0 



71 



(111.38) 



we get (note that Y results to be a real number) 

oo oo 



W a (X,fi,v) 



1 



y/if\r 



EE 

ra=0 m=0 



w.ml 



r 



m 



x a 



V2\i 



a 



y/2\i 



H n (Y) H m (Y) 



(111.39) 



This in turn must be equal to a series expansion in w nm 19 



ILK 



1 1 E rr-r w »(^/ i ' y ) 



n,m=0 



ynlrnT 



(111.40) 



so that we have 
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1 1 / * \ n / \ m 

w nm {X,fi,u) = -^--e- Y2 ^===[ r --) i-r-r) H n (Y)H m (Y). (111.41) 
V^l r l y/n\m\2 n+m \\r\) \\r\ ) 

In conclusion the marginal distribution for the n-th excited state results to be 

w n {X,p,v) = w nn {X,ii,v) = -L-e~ Y2 ^-H 2 n (Y) . (111.42) 

V^M n\2 n 

Let us compare our results with simpler cases. Assuming u(t) = u = 1 and / = 
we must recover the results of ref. [jn|. We have 



e(t) = e u , /3{t) = , (111.43) 

consequently 

Y = n = / ^ 2 ( IIL44 ) 

\r\ ^ n z + v l 



and we recover the result obtained for the simple harmonic oscillator in [19| 



To give an idea of what happens if u(t) ^ const, we show below some plots of w n (x, t) 
for different values of the parameters. In this respect it is necessary to solve Eq. ( [[1.12]) 
for e(t). In the parametric resonance case (in dimensionless units) 

9 . . 1 + k cos2t , _ * — 

co 2 (t) = — , k « 1 (HI.45) 



a good approximation for e(t) is given by |34| (see also pSj): 

e(t) = cosh ^e" - i sinh ^e - * (111.46) 

Thus, in this case, w n (x, t) is completely determined. Fig.|I] shows the plot of wo(x,t) 
when fi = 1, z/ = and = 0.01. Since in this case H (Y) = 1 we obtain the 
typical Gaussian centered around the origin of the axes, modulated in time by the factor 
l/|r| = (1 — sinh y sin t) -1 / 2 . When n increases we recognize the appearance of the 
typical structure of minima and maxima due to the zeroes of the Hermite polynomials 
Fig.f|. Quite interesting is the behaviour of w n (x,t) with respect to \i and u, which 
can be seen in Fig.|] and Fig.|| (see also Fig.|| and Fig.|[). To study the dependence of 
w on n and v we must remember that the information contained in the MDF is over- 
complete in the sense that we can choose different "tomography schemes" [0. This is 
the counterpart of the different representations (coordinate, momentum etc.) existing 
in the usual formulation of quantum mechanics. In the following we use the optical 
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tomography scheme p2lp3j: /i 2 



1. In Fig. |5| we show wq as a function of x and 



fi, we note that the maximum of probability shifts from x > to x < when fi goes 
from to 1 corresponding to the two extreme cases X = p and X = q respectively 
Fig. |^ shows the change of to at a fixed point i as a function of t and /i. We observe a 
typical oscillation behaviour at fixed \i for varying time and at a fixed instant we note 
the change of w with respect to fi. 



IV. CONCLUDING REMARKS 



In this paper we studied the driven harmonic oscillator in the framework of the 
probability representation of quantum mechanics. By means of the time-dependent 
invariants we determine the classical propagator, II, of the evolution equation of the 
marginal distribution function relative to the potential considered. In this way we were 
able to reconstruct the quantum propagator for the density matrix and, up to a phase 
factor, the quantum propagator for the wave function. We recover well known limit 
cases [[BJ] and [R3[]. 

We compute the marginal distribution function, w a , for coherent states. We obtain it 
first in the framework of the probability representation, then using the usual techniques 
of quantum mechanics. The time dependence of w a is achieved both by means of the 
classical propagator and, by directly solving a time dependent equation which encodes 
the notion of coherence. 

Starting from w a we compute the marginal distribution eigenf unctions in the energy 
eigenstate basis and study its behaviour in a particular case (the parametric resonance). 

Our results drive us to the conclusions that the quantum description of the forced 
parametric oscillator can be given in a selfconsistent approach, which is alternative to the 
Schrodinger picture while closer to the classical description. This is a further evidence for 
the possibility of formulating quantum mechanics by means of the MDF associated to a 
random variable X, avoiding complex wave functions and the density matrix formalism. 
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FIG. 1. wo(x,t) for /x = 1 and v 
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FIG. 3. wq(x, t) for fj, = and u = 1 
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